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ABSTRACT 

Spitzer and Caltech Submillimeter Observatory (CSO) images and spec- 
trophotometry of e Eridani at wavelengths from 3.5 to 350 /im reveal new details 
of its bright debris disk. The 350 /im map confirms the presence of a ring at 
r = 11-28" (35-90 AU) observed previously at longer sub-mm wavelengths. The 
Spitzer mid- and far-IR images do not show the ring, but rather a featureless disk 
extending from within a few arcsec of the star across the ring to r ~ 34" (110 
AU). 

The spectral energy distribution (SED) of the debris system implies a conplex 
structure. A model constrained by the surface brightness profiles and the SED 
indicates that the sub-mm ring emission is primarily from large (a ~ 135 /im) 
grains, with smaller (a ~ 15 /im) grains also present in and beyond the ring. 
The Spitzer IRS and MIPS SED-mode spectrophotometry data clearly show the 
presence of spatially compact excess emission at A > 15 /im that requires the 
presence of two additional narrow belts of dust within the sub-mm ring's central 
void. The innermost belt at r ~ 3 AU is composed of silicate dust. 

A simple dynamical model suggests that dust produced collisionally by a 
population of about 11 M e of planetesimals in the sub-mm ring could be the 
source of the emission from both in and beyond the sub-mm ring. Maintaining 
the inner belts and the inner edge to the sub-mm ring may require the presence 
of three planets in this system including the candidate radial velocity object. 

Subject headings: infrared: stars - circumstellar matter - planetary systems - 
stars: individual (epsilon Eridani) 

1. Introduction 

The very nearby (3.22 pc; HIP 16537) K2 V star e Eridani (HR 1084, HD 22049, GJ 144) 
was found by IRAS to have a bright semi-resolved far-infrared excess due to circumstellar 
dust (Gillett 1984; Aumann 1985). The system shows > 1 Jy of excess at 60-200 /im (IRAS 
03305-0937; Walker & Heinrichsen 2000 and references therein), and a moderate-significance 
excess at 25 /im. 

VLTI angular diameter measurements plus evolutionary models imply an age for e Eri 
of 850 Myr (Di Folco et al. 2004), consistent with previous determinations based on Ca II 
activity (Henry et al. 1996), Li abundances (Song et al. 2000), and gyro chronology (Barnes 
2007). For e Eri, as well as the other IRAS-discovered prototype debris disk systems a 
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Lyr (Vega), a PsA (Fomalhaut), and (3 Pic, (Gillett 1984; 1986) later dubbed the "Fab 
Four", estimated time scales for dust destruction by Poynting- Robertson (P-R) drag and 
mutual collisions are significantly less than the respective system ages. Thus, the observed 
dust is not primordial but must be "second-generation" material released relatively recently, 
by collisions or other activity, from larger parent bodies. These debris disks are therefore 
understood to represent late stages, or remnants, of the planet formation process. Many 
more debris disks have been discovered since the prototypes (see reviews by Backman & 
Paresce 1993 (Protostars and Planets III); Lagrange et al. 2000 (Protostars and Planets IV); 
Meyer et al. 2007 (Protostars and Planets V)). e Eri is especially interesting for its proximity, 
its age of nearly 1 Gyr, and its K main sequence spectral type for which relatively few debris 
disks are presently known (Trilling et al. 2008). 

The e Eri system was mapped at 850 and 450 /xm by Greaves et al. (1998; 2005), 
revealing a nearly circular ring extending more than 25" (80 AU) from the star, with peak 
surface brightness at about 18" (~ 60 AU) and a central cavity with radius about 10" (~ 
30 AU). Schiitz et al. (2004) detected the ring at 1.2 mm and confirmed its spatial extent. 
Marengo et al. (2006) made a deep search at wavelengths of 3.5 to 7.9 for planetary- 
mass companions between 14" (45 AU) and 2.9' (560 AU) from the star with Spitzer's IRAC 
camera, and also set upper limits to extended emission. Profntt et al. (2004) used HST/STIS 
observations to determine an upper limit to optical scattered light from the region of the 
sub-mm ring. Di Folco et al. (2007) set an upper limit on near-IR (K' band) scattered light 
brightness within 1" (3 AU) of e Eri via CHARA interferometric measurements. 

Dent et al. (2000), Li, Lunine & Bendo (2003), and Sheret, Dent & Wyatt (2004) 
calculated disk models constrained by the surface brightness distribution in the sub-mm 
images and pre-Spitzer SED data, examining alternatives of porous or solid grain structures 
and silicate or silicate-plus-ice compositions. Bright spots in the sub-mm ring have been 
interpreted as dust density enhancements caused by a planet orbiting in or near the ring 
(Liou & Zook 1999; Ozernoy et al. 2000; Quillen & Thorndike 2002; Deller & Maddison 
2004; Greaves et al. 2005). 

Tentative detection of a Jovian-mass planet orbiting e Eri with a semi-major axis of 3.4 
AU, well inside the sub-mm ring, has been made via radial velocity measurements (Campbell 
et al. 1988; Hatzes et al. 2000), possibly confirmed astrometrically (Benedict et al. 2006). 
Das & Backman (1992) suggested that such a planet might influence the location of the 
warm dust detected by IRAS. Detailed models of the gravitational effect of a planet close to 
the star on possible nearby dust distributions were calculated by Moran, Kuchner & Holman 
(2004). 

The present paper is part of the Spitzer all-instrument GTO "Fab Four" program making 
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comprehensive studies of the prototype planetary debris disks, to put our solar system into 
evolutionary context and to search for signs of planets (Fomalhaut: Stapelfeldt et al. 2004; 
Vega: Su et al. 2005; j3 Pic: Chen et al. 2007). Spitzer obtained the first far-IR images 
of e Eri, which yield broadband photometry and disk surface brightness profiles as well as 
high-resolution comparisons with prior sub-mm maps. Those Spitzer images are combined 
here with medium-resolution IR spectroscopy plus a 350 /im map from the CSO SHARC 
II camera. Together, these data provide a unique opportunity to study the nearest known 
exo-Kuiper belt and exo-zodiacal dust system in detail. 

Section 2 of this paper describes the Spitzer and CSO observations of e Eri and the data 
reduction. Section 3 addresses detailed analyses of the data including first-order corrections 
of image and spectrophotometric data plus determination of the photosphere SED and IR 
excess. Section 4 describes derivation of a model of grain spatial and size distributions 
fitted to the disk's multi-A surface brightness distributions and over- all SED. Section 5 
discusses implications of the model results, compares them with other observations and SED- 
based models of the e Eri disk, and presents results of a simple dynamical model calculation 
connecting e Eri's disk properties with the history of our solar system's Kuiper belt. 



2. Observations and Data Reduction 

Spitzer observations of e Eri were taken in Guaranteed Time Observing Program 90 
(P.I. Michael Werner). Data obtained with various instruments and facilities were used 
as follows: (1) Spitzer IRAC image photometry measures the photospheric SED from 3.6 
to 7.9 /im, wavelengths that show no detectable contribution from the debris disk. (2) 
Spitzer MIPS data at 24, 70, and 160 /im provide mid- and far-IR images of the disk. (3) 
Spitzer MIPS SED-mode data give source flux densities from 55 to 90 /im with low spectral 
resolution, plus source size in one dimension versus wavelength. (4) Spitzer IRS spectra 
link the mid-IR broadband photometric points, determine the shortest wavelength of excess 
(highest abundant dust temperature, minimum dust orbit radius), and allow a search for 
mineralogical features. (5) The CSO 350 /im map can be compared with the Spitzer far-IR 
images as well as with prior 450 and 850 /im maps. Surface brightness profiles from the 24, 
70, 160 and 350 /im data together strongly constrain the disk model. (6) The Spitzer MIPS 
images and CSO map are also integrated to yield broadband photometric points for the 
entire disk, supplemented by IRAS catalog photometry using re-calculated color corrections, 
plus previously published ISO photometry. 



- 5 - 



2.1. Spitzer IRAC data 

Images at 3.6 to 7.9 of e Eri and its surroundings were made with Spitzer's InfraRed 
Array Camera (IRAC; Fazio et al. 2004) in 2004 January and February. Those observations 
were designed to look for wide-separation low-mass companions, but also provided accurate 
photometry of the star via Point Spread Function (PSF) fitting of non-saturated portions of 
e Eri and standard star images. Details of the observations and data reduction procedures 
are described by Marengo et al. (2006). The IRAC data were reprocessed for the present 
paper to take advantage of an updated instrument calibration plus availability of additional 
reference star observations to generate a new PSF. 

Each of the two IRAC observations consisted of a 36-position dither of the source on 
the array, for a total of 2 x 3380 seconds integration time. The new PSF reference stars were 
e Indi (K4.5 V), Vega, Fomalhaut and Sirius. (The debris disks around Vega and Fomalhaut 
are not detected at IRAC wavelengths.) The total field of view had a width of 5.78' due to the 
dither, slightly larger than the IRAC frame width of 5.21'. Data reduction and calibration 
were done using the Spitzer Science Center (SSC) pipeline version S15 combined with the 
post-BCD IRACproc software (Schuster, Marengo & Patten 2006). 

2.2. Spitzer MIPS data 

2.2.1. MIPS 24, 70 and 160 fim Images 

Observations of e Eri were made using all three imaging channels of the Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004). Images at 24 /im were made on 
2004 January 29 using two different observing strategies. (1) The standard 24 fxm MIPS 
photometry dither pattern was executed at four cluster target positions in a square pattern 
8.29" (3.25 pixels) wide. This procedure yielded an exposure time of 4x(16x3) seconds for 
the region within 2.5' of the star, with a high degree of sub-pixel sampling to facilitate PSF 
subtraction. Observations of the K4 III star HD 217382 made on 2003 November 19 with 
the identical dither/cluster target pattern were used as a PSF reference. (2) To probe a 
wider field of view for possible companions, another e Eri 24 /im dataset was taken using a 
3x3 raster map. The resulting mosaic covers a 0.25° xO. 25° region centered on the star to 
an exposure depth of 2 x (16x3) seconds. 

Two types of MIPS 70 /im image observations were also taken on 2004 January 29. (1) 
With the default pixel scale (9.96"/pixel), eight cycles of the large- field photometry dither 
pattern were used, providing a total exposure time of 151 seconds per pixel. (2) With the 
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fine pixel scale (5.2"/pixel), one cycle of the small field photometry dither pattern was made 
at each of 12 cluster target positions, arranged in a rectangular grid with 16.22" (3.25-pixel) 
spacing, for an effective exposure time of 12 x (10x3) seconds per pixel. On 2005 April 4 a 
deeper fine scale image was obtained with an effective exposure time of 12 x( 10x10) seconds 
per pixel. 

Two MIPS 160 /im image datasets were obtained. (1) The first set was taken on 2004 
January 29, using 2 cycles of small-field photometry at each of 9 cluster target positions 
spaced on a 36" rectangular grid, and 3-second exposures, for an exposure time of 9x27 
seconds per pixel. (2) Another set was obtained on 2006 Feb 19 using the same observing 
strategy except that 10-second exposures were used, for an exposure time of 9x90 seconds. 

The MIPS 160 camera suffers from a spectral leak by which stray light in the 
wavelength range 1-1.6 jum produces a false image partially overlapping the actual 160 /xm 
source image. This ghost image is several times brighter than the real disk emission in the 
case of e Eri, and must be subtracted accurately to enable study of the circumstellar emission. 
The spectral leak is field-dependent, and therefore can only be properly subtracted using a 
reference star that was observed with the same dither parameters as the science target. Pre- 
launch observation planning did not anticipate this requirement, so no appropriate reference 
star observation was available initially to leak-subtract the e Eri dataset. 

Good leak subtraction was achieved using 160 fxm images taken on 2004 August 6 of the 
debris disk star r Ceti (type G8) with the same dither pattern as for e Eri. Only the spectral 
leak, and not the photosphere or disk excess, is detected at 160 //m from r Ceti (Stapelfeldt 
et al., in preparation). The subtraction procedure involved empirically determining the max- 
imum normalization factor for the leak reference source such that its subtraction from the 
science target did not produce noticeable residuals below the background level. The subtrac- 
tions were performed using mosaics in detector coordinates of all the dithered observations. 
Observations of Achernar (a Eri) were also obtained for e Eri leak subtraction reference, 
on 2005 August 30. The Achernar data did not support a good leak subtraction for e Eri, 
probably because of the large color difference between the two stars (type B3 versus K2). 

The MIPS instrument team's Data Analysis Tool (DAT; Gordon et al. 2005) was used for 
basic reduction (dark subtraction, flat fielding/illumination correction) of the SSC Sll data 
products. The 24 /zm images were processed to remove a vertical "jailbar" pattern along the 
detector columns that appears during observations of bright sources. The known transient 
behaviors associated with the MIPS 70 //m array were removed from the default-scale data 
by median column subtraction and time-filtering the images with an excluded region of 15 
pixels (~ 150") centered on the source (for details see Gordon et al. 2007). For the 70 /zm 
fine-scale data, off-source exposures were subtracted from preceding on-source exposures to 
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remove the time- dependent background column offsets. Flux calibration factors of 0.0454 
MJy sr" 1 (24 /mi), 702 MJy sr" 1 (70 /iin default-scale), 2894 MJy sr" 1 (70 /mi fine-scale), 
and 44.7 MJy sr" 1 (160 /mi) per MIPS raw data unit were applied to the MIPS image data 
(Engelbracht et al. 2007; Gordon et al. 2007; Stansberry et al. 2007). The processed data 
were then combined using World Coordinate System (WCS) information to produce final 
mosaics with pixels half the size of the physical pixel scale. 

2.2.2. MIPS SED-mode Data 

MIPS Spectral Energy Distribution (SED-mode) data were obtained on 2005 September 
4 using 10 cycles of 10 second exposures with a +3' chopper throw. The total integration 
time was 629 seconds. Basic reduction of the data to produce calibrated versions of the 
individual exposures was done using version 3.06 of the MIPS DAT code. MIPS SED- 
mode observations provide chopped pairs of images. The off-source images were pair-wise 
subtracted from the on-source images, and the results mosaicked to produce an image of the 
55-90 /mi spectrum of e Eri. The mosaic was generated using 4.9" pixels, yielding an image 
65 pixels across in the spectral dimension by 44 pixels in the spatial dimension. This image 
was then boxcar-smoothed over 5 adjacent rows, resulting in spectral resolution R = 15-20. 

2.3. Spitzer IRS data 

Mid-infrared spectra of e Eri were obtained in 2004 February, 2004 August, and 2005 
January with the Spitzer Infrared Spectrograph (IRS; Houck et al. 2004). Using the low- 
spectral resolution (R ~ 60-120) spectrographs SL1 (7.5-14 /mi), LL2 (14-21 /mi) and LL1 
(20-38 /mi), observations were made with five slit positions separated by half a slit width 
(1.8") in SL1 and ten slit positions separated by half a slit width (4.8") in LL2 and LL1. 
The star was centered in the scan pattern in each case. Staring-mode observations were also 
obtained using the low-resolution SL2 (5.3-7.5 /mi) and high resolution (R ~ 600) SH (10-20 
/mi) and LH (19-37 /mi) modules. Those observations were at two nod positions, with the 
star displaced from the center of the slit by a third of the slit length. Total integration times 
per slit position or nod were 84 seconds in the SL1 spectral map, 70 seconds each in the LL2 
and LL1 spectral maps, and 54 seconds, 72 seconds, and 144 seconds respectively for the 
SL2, SH and LH observations. Each set of observations was preceded by a high-accuracy 
pointing peak- up observation of the neighboring star SAO 130582 using the IRS Blue (15 
/mi) camera. 
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Data reduction began with products of the Spitzer Science Center's Sll IRS data 
pipeline. Further processing involved removal, by interpolation in the spectral direction, of 
permanently bad and "rogue" pixels identified in the IRS dark-current data for all campaigns 
up to and including the one in which each e Eri spectrum was taken. The IRS instrument 
team's SMART software was then used to extract 1-D spectra from the 2-D long-slit (SL, 
LL) or echelle (SH, LH) images. Sky emission in the SL and LL observations was measured 
by a linear fit to the signals at points along the slit farther than 40" from the star, and 
subtracted as part of the extraction. The e Eri system appeared to be spatially unresolved 
in the SL observations, so an extraction window was used that was fitted to the instrumen- 
tal point-spread function at the longest wavelengths of each low-resolution module, 4 pixels 
wide in each case. A full-slit extraction was performed for the high-spectral-resolution (SH 
and LH) observations. The sky emission measured in the low-spectral-resolution observa- 
tions was small enough compared with the flux received from e Eri and its disk that sky 
subtraction in the SH and LH observations could be neglected. 

Reference SL and LL observations of a Lac (Al V) plus SH and LH observations of 
£ Dra (Kl III) were obtained in the same manner as the e Eri observations. The spectra 
of e Eri were calibrated via division by spectra of the respective reference star taken at the 
corresponding nod position, multiplication by a synthetic template spectrum of the reference 
star (Cohen 2004, private communication), and tben nod- averaging. Reduced spectra from 
each spectrograph module were then combined to produce complete 5.3-38 /im spectra. 
Small spectrograph module-to-module differences at wavelengths of overlap were assumed 
to be due to minor telescope pointing-induced throughput losses, and removed by scaling 
up the spectra from smaller-slit modules to match those from the larger-slit modules in the 
overlap range. 

Linearity and saturation in the IRS detector arrays is corrected in the data pipeline, 
but e Eri and its disk are bright enough that linearity and saturation effects in the low- 
spectral-resolution modules are large and the corrections uncertain. The main impact of 
these effects is that the LL spectra had significantly lower effective signal-to-noise ratio and 
spectrophotometric accuracy than the SH and LH observations due to residual inaccuracy 
of the LL saturation correction. The final spectrum used in this paper was assembled from 
SL, SH and LH data; LL data will not be discussed further. The uncertainty in flux density 
is 5% below about 10 Jy, increasing to about 10% for flux densities as high as 20 Jy. 
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2.4. CSO Data 

Observations of e Eri at 350 /im were made during parts of three nights in 2003 January 
and two nights in 2005 December using the SHARC II 12x32 bolometer array (Dowell et 
al. 2003) at the Caltech Submillimeter Observatory (CSO). The weather conditions were 
excellent, with 225 GHz (1.3 mm) zenith opacity ranging from 0.028 to 0.048. The total 
integration time on the source was 11 hours. Data were obtained by rapid scanning (without 
chopping) in Lissajous patterns with periods of 14-20 seconds and amplitudes of ±15" to 
±40" in each of two axes. 

Each SHARC II detector spans 4.6", and the diffraction-limited beam size is 9" (FWHM), 
with a PSF determined by coadding images from a variety of calibration sources interspersed 
with the e Eri observations. Data reduction used the sharcsolve software package that iter- 
atively decomposes fixed source emission, time variable atmospheric emission, and detector 
offsets. Best results for subtraction of the temporally and spatially variable telluric emission 
were achieved by assuming that the average emission outside a 35" radius from the star 
was zero, at the cost of filtering out possible more widely extended emission. Calibration 
was based on the integrated map signals within a 35" radius, following correction for atmo- 
spheric absorption, from sub-mm standards OH231.8 (18 Jy), CRL618 (19 Jy), o Ceti (2.3 
Jy), Ceres (13 Jy in 2003 January), and Juno (8.0 Jy in 2005 December). The uncertainly 
in the absolute flux scale is estimated to be 25%, and the RMS noise level in the reduced 
350 /im image is ~ 3 mJy per smoothed beam. 



3. Results 
3.1. Images and Spatial Structure 

3.1.1. MIPS Far-IR Images 

The top row of Figure 1 displays the MIPS 24, 70 and 160 /im images after standard 
Spitzer pipeline processing. The bottom row of Figure 1 shows the output from a further 
level of processing. Specifically: (1) Figure Id is the result of subtracting a scaled PSF from 
the 24 /jm image in la, (2) Figure le shows a 70 /im fine pixel-scale image for comparison 
with the default pixel-scale image in lb, and (3) Figure If is the result of compensating for 
the 160 /im filter's short-wavelength leak apparent in lc. Radial profiles of the processed 
images are presented and discussed in Section 4.2.1. 

The 24 /im PSF-subtracted data show an image core size of 5.6" x 5.4", consistent with 
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the values of 5.52" x5. 48" for an unresolved blue (photospheric) source determined by Su et al. 
(2008). The stripes in the processed data are PSF subtraction residuals, mainly noise arising 
from different sampling of the source and PSF reference by the pixel grid. Examination of 
the 24 /im PSF wings shows evidence for slight excess surface brightness at the position of 
the first Airy dark ring (r ~ 7") but with only ~ la significance. 

The far-IR source is resolved at 70 and 160 /im. At 70 /im (fine pixel-scale image) 
the observed source has a compact core with FWHM of 24.2" x 22.3" elongated approxi- 
mately north-south, superimposed on a relatively flat pedestal of emission roughly 1' in 
diameter. The core 70 /im emission lies mostly interior to the sub-mm ring, and is clearly 
extended relative to the size of unresolved point sources, which is 15.3" x 14.4" (Su et al. 
2008). The leak-subtracted 160 /im source can be approximated by Gaussian profiles with 
FWHM 62"x52", also elongated north-south. This distribution is also extended relative to 
a point source (39" diffraction FWHM). The stellar point source contributes about 10% of 
the total system flux density at 70 /im and 5% at 160 /im. 

The 70 and 160 /im images contain a substantial surprise. Although the far-IR source 
FWHM is about the size of the sub-mm ring, the far-IR images do not show the ring- 
with-deep-hole aspect observed at longer wavelengths. Instead, the far-IR surface brightness 
distribution is continuous across the central region and the sub-mm ring, extending well 
beyond the ring, with no significant enhancement at the position of the ring. The far-IR 
source is close to circularly symmetric, with possible slight extension north-south as in the 
350 and 450 /im source maps, but unlike the 850 /im maps that are extended east-west or 
northeast-southwest (Greaves et al. 2005). 

To test for far-IR detection of surface brightness features seen in the sub-mm ring, an 
azimuthal average of the 70 /im fine-scale image was subtracted from the image itself. Figure 
2a shows the result of this analysis. The subtracted image is smooth to < 20% of the surface 
brightness in the ring, smaller than the amplitude of features in prior sub-mm maps. The 
70 /im image therefore does not show the azimuthal surface brightness variations observed 
at longer wavelengths. A similar result, albeit with lower signal-to-noise ratio, is obtained 
with the 160 /im image (Figure 2b). 

3.1.2. MIPS SED-mode Source Size 

Figure 3a shows the MIPS SED-mode source FWHM (along the slit) versus wavelength 
from Gaussian fits to the spatial profiles of the signal at each row of the mosaic image. The 
FWHM are 21, 23.5, and 38" at 55, 70 and 90 /im, respectively, consistent with the 70 and 160 
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/im image source sizes. Analysis of spectrophotometric data from the MIPS-SED instrument, 
presented in Figure 3b, is described in Sections 3.3.1 and 4.3.2. The small-scale variations in 
FWHM and in flux density versus wavelength relative to simple monotonic trends, as well 
as the steep increase in FWHM longward of 87 /im, are not real but represent systematic 
uncertainties in the reduction procedures, which dominate photometric uncertainties. 

3.1.3. CSO Sub-mm Map 

Figure 4 displays the CSO SHARC II 350 /im map. The sub-mm ring that has been 
observed previously at 450 and 850 /im is apparent here, with approximately the same spatial 
scale. The 350 /im image has flux consistent with zero excess in the central beam (FWHM 
= 29 AU) after subtraction of a model 8 mJy point source corresponding to the stellar 
photosphere emission, representing 2% of the total flux in the map (Table 1). The 350 jum 
ring may show north-south extension, as in the Spitzer 70 and 160 /im images and the JCMT 
450 /im map, but differing from the 850 /im morphology. 

The average value and dispersion of the surface brightness around the ring at 350 /im is 
S u = 0.15 ± 0.04 mJy per square arcsec. The brightest 350 /im spot, in the north to north- 
east sector of the ring, has S^ ~ 0.24 mJy per square arcsec, about 2a above the average. 
The location and shape of the 350 /im ring are confirmed by the 2003 and 2005 data sets 
analyzed independently, but the azimuthal substructure is not. 

3.2. Broadband Photometry 

3.2.1. Stellar Photosphere 

Precise determination of e Eri's photospheric SED was necessary to measure the excess 
emission from the disk, especially in the mid-IR where the system SED is only slightly above 
the stellar continuum. The star has m^ = +3.74, B-V = +0.88 (Mendoza et al. 1978) and 
V-K of +2.11 (Koornneef 1983a). Average main sequence B-V colors are +0.82 and +0.92 
for K0 and K2 temperature classes, respectively (Johnson 1966). Average V-K colors are 
+2.00 and +2.25 for Kl and K2, respectively (Koornneef 1983b). These data imply that e 
Eri is a little bluer than the average K2 V star, with interpolated temperature class between 
Kl and K2. Flower (1996) gives T e// = 5090 K for B-V = 0.88. 

Flux densities from e Eri at 3.6, 4.5, 5.7 and 7.9 /im determined by PSF-fitting of Spitzer 
IRAC data (Section 2.1) have ratios nicely consistent with expectations for a photosphere 
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of this star's temperature and luminosity. The photosphere model SED used in this paper 
is a Kurucz model for T = 5000 K (nearest match to 5090 K) and log g = 4.5, normalized 
to the four IRAC flux densities via least-squares minimization. The RMS of the fit to those 
points is 1.5%. The result is consistent with ground-based photometry at 2.2 fim and shorter 
wavelengths to better than 2% precision. Figure 5 displays the IRAC data points and fitted 
photosphere in relation to other broadband measurements discussed below. Table 1 presents 
extrapolated photospheric flux densities which can be regarded as having 2% uncertainty. 

3. 2. 2. Photometry from Image Data 

Since the disk is not resolved at 24 fim, the total emission from the disk at that wave- 
length was obtained using both aperture photometry and PSF-fitting photometry optimized 
for point sources; for details see Su et al. (2006). A total of 2.04 ± 0.04 Jy was estimated 
for the star and disk combined, which is 18% above the expected photosphere. No color 
correction was applied at 24 fim because the observed SED slope in that wavelength region 
implies a MIPS filter correction factor of ~ 1. 

The total flux density of the system at 70 and 160 fim was integrated over a circular 
aperture with radius 64" centered at the star. The total disk emission in the 70 /im and 
160 fim bands is respectively 1.50 and 0.93 Jy after application of color correction factors of 
0.893 and 0.971 appropriate for a 55 K blackbody. 

IRAS photometry of e Eri spans an important range of wavelengths, from the red end 
of the photospheric SED to beyond the peak of the excess. The Spitzer IRS spectrum, MIPS 
images and MIPS SED-mode data allow a refined estimate of the continuum slope across the 
IRAS 25 and 60 fim bands and hence more accurate color corrections, which were applied 
here. 

Figure 5 and columns 1 & 2 of Table 1 present the merged system (star+disk) broadband 
SED including Spitzer MIPS photometry from image data, integrated flux density from the 
CSO 350 fim map, newly color-corrected IRAS catalog flux densities, ISO 200 fim flux 
density, and previously published integrated flux densities at 450 and 850 fim from Greaves 
et al. (2005). The peak of the SED is between A ~ 70 and 100 fim. 

ISO photometry of e Eri (Walker & Heinrichsen 2000) yielded a flux density of 1880 
mJy at 200 fim (cf. predicted photosphere value of 25 mJy), far above the 160 fim point. 
The Spitzer 160 fim image, as well as the 450 and 850 fim sub-mm maps, show significant 
background emission extended beyond 40" from the star that would have been included in 
the large long-wavelength effective beam sizes of both ISO (200 fim photometric aperture 
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diameter 90") and IRAS (100 /im detector footprint ~ 180" x 300"). Numerous objects with 
near-IR colors of galaxies that would be strong sources at A > 100 fim are detected in 
Spitzer IRAC images of the e Eri field (Marengo et al. 2006). The presence of this complex 
background emission probably accounts for the high flux density measured by ISO; the 
present observations indicate that much of the ISO flux must come from beyond the sub-mm 
ring. The IRAS 100 /im and ISO 200 /im photometry (Figure 5) therefore provide upper 
limits to the system flux densities at those wavelengths. A less confused view of the source 
at 100-200 /im should be possible eventually with deep Herschel PACS and SOFIA HAWC 
observations. 



3.3. Spectrophotometric Data 

3.3.1. MIPS SED-mode 

With e Eri expected to appear as an extended source along the MIPS SED slit, simple 
aperture photometry measurements cannot be made without admitting significant back- 
ground noise. To address that problem, a Gaussian was fit to the spatial profile of the signal 
along each row of the SED mosaic. The integrated product of the Gaussian width and height 
gives the source response at each wavelength. The resulting e Eri spectrum in instrumental 
units was converted to physical units by reducing and extracting an archival spectrum of 
the calibration star Canopus (a Car, G8 III) in the same way. The Canopus spectrum was 
normalized to give a flux density of 2.15 Jy at 70 fxm (the star's 3.11 Jy source flux density 
multiplied by a 69% slit throughput factor appropriate for a point source). The result was 
then corrected for the spectral response function by assuming Canopus has a Rayleigh- Jeans 
spectral slope in the far-IR. The flux normalization and spectral response function derived 
from Canopus were then applied to the e Eri spectrophotometry. 

Figure 3b presents the resulting MIPS SED-mode spectrum from 55 to 90 /im, giving 
the detected flux density transmitted through the MIPS SED slit independent of assump- 
tions about the true spatial extent of the source. The slope of the spectrum, declining 
to longer wavelengths, is due to the intrinsic shape of the star+disk SED combined with 
wavelength-dependent slit losses as the size of the diffracted and extended emission region 
varies relative to the 20"-wide MIPS SED entrance aperture. Linking the MIPS SED-mode 
spectrophotometric data to a disk model is described in Section 4.3.2. 
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3.3.2. IRS 

Spitzer IRS spectrophotometric data for the disk plus stellar photosphere are included 
in Figure 5. The photosphere model described in Section 3.2.1 is overplotted for comparison. 
A constant offset of 173 mJy was subtracted from the IRS data, determined by making the 
average excess above the photosphere model from 5 to 12 /iin equal zero. This is consistent 
with the observed lack of significant excess emission in bradband photometry at 5.7, 7.9 
and 12 /im (Table 1). That offset represents only 1-3% of the total source flux in that 
wavelength range, well within the IRS calibration uncertainty. The resulting rescaled IRS 
mid-IR SED is in close agreement with the MIPS 24 /im image photometry as well as with 
the lower-precision IRAS 25 /im broadband point. 

The Figure 5 panel inset shows the spectrum from 10 to 30 /im after subtraction of the 
photosphere. Significant departure of the e Eri SED from the photosphere begins at about 
15 /im. The mid-IR SED is remarkable in having a very steep rise from 15 to 20 /im, a 
"plateau" of almost constant flux density from about 20 to 30 /im, and then the beginning 
of another steep rise beyond 30 /im. 

4. Modeling 

4.1. General Considerations 

The data presented above provide new information on the e Eri circumstellar dust spatial 
distribution, new spectral information regarding warm dust components, and separation of 
disk emission from confusing background sources. A model of particle spatial distributions 
and physical properties can now be produced that will allow improved understanding of 
dust lifetimes and dynamics, which in turn allows inferences about locations of planetary 
companions, and about the evolution of this system for comparison with our solar system. 

The modeling process involved: (1) derivation of intrinsic 70, 160 and 350 /im 1-D 
radial surface brightness profiles that, when convolved with the respective PSFs, match the 
observed radial profiles; (2) calculation of a physical model for the outer disk with specific 
grain spatial distributions, size distributions and compositions that yields the intrinsic 70, 
160, and 350 /im surface brightnesses inferred from the first step; (3) derivation of another 
physical model for warm material detected by IRAS and Spitzer IRS but not resolved in 
the imaging data, then testing that the predicted long-wavelength emission from the inner 
material would be consistent with the 70, 160 and 350 /im radial profiles; and finally, 4) 
checking that the over-all physical model, including both resolved and unresolved material, 
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correctly predicts the observed spectral energy distribution of IR excess from 15 to 850 /im 
and also the MIPS SED-mode narrow-slit flux densities plus 1-D source sizes or limits. As 
mentioned in Section 3.3.1, the MIPS SED fluxes need to be corrected for slit losses, and 
that can only be done by estimating the corrections iteratively using synthetic images derived 
from the model (Section 4.3.2). 

The 350 jum map shows a deficit of sub-mm emission at r < 35 AU relative to 35-90 
AU, similar to previous sub-mm maps. In contrast, the 24, 70 and 160 /im image profiles and 
photometry plus mid-IR spectrophotometry and broadband photometry, as well as the orig- 
inal IRAS detection, together indicate the presence of significant amounts of warm material 
lying at r < 35 AU, within the central void in the sub-mm ring. 

The surface brightness in the high-S/N 70 /im images is continuous and monotonic out 
to r > 45 arcsec (140 AU), far outside the sub-mm ring. However, trial PSF convolutions 
of model intrinsic brightness distributions revealed that the outer boundary of the 70 jum 
emission is actually located at r ~ 34" (110 AU), with apparently further extended emission 
representing the PSF wings. In comparison, nearly all of the observed sub-mm surface 
brightness is contained in the ring between about 35 and 90 AU (11-28") (Greaves et al. 
2005; also, see Schiitz et al. 2004 regarding the source diameter at 1.2 mm). The region 
outside the sub-mm ring, beyond r = 28" (90 AU), is here designated the "halo". 

An outer limit to the 160 /im emission is more difficult to determine due to background 
source confusion discussed above, but the data and preliminary modeling indicate that the 
over-all spatial distribution of material observed at 160 /im much more closely resembles the 
70 /im than the 350 /im emission. Significant emission at 350 /im does not extend beyond 
28" (90 AU). 

4.2. Outer Fully Resolved Disk 

4-2.1. Radial Surface Brightness Profiles 

The 70 /iin images have small ellipticity and also are azimuthally smooth (Figure 2), 
justifying uses of 1-D average radial profiles centered on the star position for detailed com- 
parison of model surface brightness distributions to observations. Based on inspection of the 
radial profiles and preliminary modeling, fiducial radii of r = 11, 28, and 34" (35, 90, and 
110 AU), representing transition points in some of the surface brightness distributions, were 
chosen as primary points for matching models to observations. 

Intrinsic surface brightness profiles at 70, 160 and 350 /im were constructed by choosing 
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surface brightness values at the three fiducial radii and interpolating between those points 
with simple power laws. Production of model PSF-convolved radial profiles for trial-and- 
error comparison with the observations involved expansion of model 1-D radial profiles into 
2-D pseudo-images, convolution of each pseudo-image with the respective wavelength's 2- 
D PSF derived for a point source with e Eri's temperature using STinyTim (Krist 2006), 
and then collapsing the convolved model 2-D map back into a 1-D profile. After each trial, 
differences between predicted and observed surface brightness profiles were used to estimate 
new values at the fiducial radii for the next trial. 

The results of modeling the 24, 70, 160 and 350 fim surface brightness profiles are 
displayed in Figure 6. The actual locations of changes in profile slopes near the fiducial radii 
are uncertain to probably 10% in radius, i.e. the modeling process had approximately that 
much spatial resolution. Note that these intrinsic surface brightness profile determinations 
were aimed merely to match the observed profiles, without reference at this step in the 
modeling process to physically motivated grain properties, which are treated in the next 
section. 



4-2.2. Physical Grain Models 

The Heidelberg DDS code (Wolf & Hillenbrand 2005) was used to estimate grain prop- 
erties and surface densities that would yield the surface brightness values determined in the 
preceding modeling step. Grains were assumed to be solid (non-porous) and spherical. DDS 
settings represented e Eri as a 5090 K blackbody with 0.33 L . Choice of grain composi- 
tion, grain minimum and maximum sizes, size distribution power law exponent, and mass 
surface densities at the three fiducial radii allowed calculation of emitted surface brightness 
for comparison with the intrinsic 70, 160 and 350 /im surface brightness values determined 
at those locations. Integrated disk flux densities and grain masses could then be calculated 
using surface densities and grain properties interpolated between the values at the fiducial 
radii. 

The large difference in appearance of the sub-mm maps versus the Spitzer 70 and 160 
/im images led to a search for model disk components such that: (1) material having high 
sub-mm and low far-IR emissivity would dominate in the sub-mm ring region, coexisting 
with (2) material having low sub-mm and high far-IR emissivity extending across the ring 
into the halo region. Many combinations of grain properties, size distributions and spatial 
distributions can match those constraints. A simple model can be built of components that 
separately fit the 70 & 160 /im images and the 350 jum map. 
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The final model for the resolved disk combines: (1) large (radii a = 100-200 jum) amor- 
phous H 2 ice grains located in the sub-mm ring, plus (2) small (a = 6-23 jum) "astronomi- 
cal" silicate grains having a broader spatial distribution that includes the ring. Experiments 
varying the value of the grain size distribution exponent indicated that the conventional 
collision equilibrium value of x — —3.5 could be used for both populations although other 
values are not excluded. Net emissivities integrated over the size distributions are approx- 
imately equivalent to single-size particles with radii of 135 /im and 15 /im for the large icy 
and small silicate grains, respectively. 

4.3. Two Inner Warm Belts 

4-3.1. Dust Properties and Locations 

Spitzer IRS, MIPS, and MIPS-SED data reveal the presence of material within the sub- 
mm ring's central void. PSF analysis of 70 and 160 /im image surface brightness profiles 
shows emission closer to the star than 11" (35 AU) but also reveals that the emission does 
not extend all the way to r = 0. The local maximum and plateau in the IR excess around 
20-30 /im indicate emission at T ~ 100 — 150 K, confirming the location of material at 
r < 35 AU, because at 35 AU, Td us t would be only ~ 40 K for blackbody grains or ~ 60 K 
for l-/im silicate grains. Quadrature subtraction of a point source size from the MIPS-SED 
1-D source extent at 55 /im (Figure 3a) yields an emission radius of very roughly 8" (25 AU). 

Modeling the properties of the inner material involved first noting that the marked 
change in SED slope at about 20 /im (Figure 5) is easily explained as a silicate mineralogical 
band (e.g. Chen et al. 2006). Trial-and-error experiments with the Heidelberg DDS code 
showed that the detectability of the band sets a rough upper limit on grain radius a < 20/2-7T 
~ 3 /im. In contrast, "generic" grains of that size without mineralogical attributes produce 
an emission SED that is significantly broader than observed, as do ice grains. 

The 15-30 /im segment of the e Eri SED can be matched approximately by emission 
from a belt of warm "astronomical" silicate grains with characteristic size ~ 3 /im, located 
about 3 AU (~ 1") from the star. This innermost warm belt is responsible for > 95% of 
the system's excess emission at 25 /im first detected by IRAS. Note that dust located 3 AU 
from e Eri would be entirely contained within the IRS slits that have a minimum width of 
4.7" (15 AU). 

The precise position of the "3 AU" belt is constrained by relations between grain size, 
radiative efficiency, and temperature. For example, if the model grain size a is set at 2.2 /im 
instead of 3 /im, the belt must move to r ~ 3.8 AU to have the correct color temperature. 
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However, grain efficiency declining steeply to longer wavelengths for such small grains would 
unacceptably "erode" the 20-30 /im SED plateau. Thus, the grain size and belt location are 
determined to perhaps 25% precision. There is no separate information regarding an upper 
limit to the belt's width, but that can be taken to be equivalent to the uncertainty in belt 
location, i.e. ~ 25%. 

The rising IRS spectrum near 35 /im (Figure 5) cannot be fit by emission from the 
material inferred at ~ 3 AU, which has an SED decreasing beyond 30 /im. The steep 
upward slope from 35 to 60 /im, characterized by a color temperature of about 55 K, plus 
the relatively narrow SED peak at 70-100 /im, locate the next lower-temperature population 
of grains also in a belt with narrow ranges of temperature and particle size. The precise 
properties of this material are especially constrained by the MIPS SED-mode 55-90 /im flux 
densities, but because of slit losses, those data must be renormalized using models of the 
source geometry, as described in the next section. 

A belt of a ~ 8/im "astronomical" silicate grains located at r ~ 20 AU (~ 6") satisfies 
the SED and source size constraints. Models placing the belt at 20 AU can be distinguished 
from alternate locations at 10 or 30 AU. Material at 20 AU could also explain the slight 
filling-in of the first Airy dark ring at 24 /im, r ~ 23 AU (Figure 6a). Silicate grains of 
this size and temperature would not show mineralogical features in the Spitzer data, so the 
composition of this belt is consistent with, but not determined to be, silicate. 

Models with continuous distributions of material extending across the sub-mm void 
produced too much total mid- and far-IR flux and/or did not match the surface brightness 
profiles and source size measurements. The warm material is clearly restricted to annular 
zones with relatively empty gaps. 

4-3.2. Modeling Details, Including MIPS SED Normalization 

Convergence on final models for the two inner warm belts involved iterative comparison 
of model SEDs from the Heidelberg DDS software with the IRS 15-35 /im and 55-90 /im 
(MIPS SED-mode) flux densities not already accounted for by the model of the resolved disk. 
Trial models yielded intrinsic surface brightness profiles at the observed wavelengths that 
were: (1) convolved with the 24 /im PSF to compare with the observed unresolved source, 
(2) convolved with a derived Spitzer 55 /im PSF and MIPS SED-mode slit mask to compare 
with the observed source 55 /im 1-D FWHM, and (3) added to the already-established 70, 
160 and 350 /im intrinsic surface brightness profiles of the disk and then PSF-convolved to 
check that the profiles were not significantly perturbed by the extra warm material. 
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The model properties of these two belts within the sub-mm void interact, in that they 
both produce significant flux density at A = 30-35 /im. For example, changing the second 
warm belt's properties toward smaller spatial radius and/or particle radii produces greater 
30-35 /im model emission, requiring the innermost belt to have greater spatial radius and/or 
particle radii so as to preserve the total mid-IR SED shape. 

The slit losses and flux normalization of the MIPS SED data were derived by generating 
2-D synthetic images of the source using model radial profiles of the disk computed at the 
MIPS SED wavelengths. This analysis was restricted to the two extremes of the MIPS SED 
wavelength range, 55 and 90 /im, using the midpoint at 70 fim as the flux normalization 
reference. 

Synthetic 2-D images at 55, 70 and 90 /im were derived by convolving azimuthal pro- 
jections of the radial disk model profiles with the STinyTIM PSF, and then adding a point 
source rescaled to the stellar flux at each wavelength. The resulting 70 /im model "image" 
closely matches the 70 /im observations (Figure 6b). A software mask with the size of the 
MIPS SED slit was then centered on the synthetic images to measured the fraction q v of the 
total flux recovered through the SED-mode slit: 58%, 47% and 38% respectively at 55, 70 
and 90 /im. The overall flux normalization factor N^ ux at 70 /im is derived by rescaling the 
SED in Figure 3b to have a 70 /im flux density identical to the MIPS 70 /im image total 
flux density (1.688 Jy, Table 1) which yields N^ x ~ 1.78. The correction for the other two 
wavelengths can be derived by factoring out the slit loss at 70 /im and multiplying by the 
slit losses at 55 or 90 /im: 



N fiux = N fiux x (1) 

This procedure produced corrected MIPS SED flux densities of 1.62 ±0.25 Jy at 55 /im and 
1.64 ±0.25 Jy at 90 /im. Uncertainties were estimated by quadrature addition of the 70 /tm 
uncertainty to the uncertainties in the uncorrected flux densities plotted in Figure 3b. Once 
the photospheric contribution is subtracted, an excess from the disk of 1.30 ± 0.25 Jy at 55 
/xm and 1.52 ± 0.25 Jy at 90 /tm is obtained. 



4.4. Model Summary 

As described in the preceding section, the over-all model of e Eri's circumstellar material 
includes: (1) particles with low far-IR emissivity and high sub-mm emissivity, consistent with 
the properties of radius a = 100-200 /im (effective a ~ 135 /im) amorphous H 2 ice grains 
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in the sub-mm ring at r = 35-90 AU, (2) particles with high far-IR emissivity and low 
sub-mm emissivity, consistent with the properties of a = 6-23 /zm (effective a ~ 15 jum) 
"astronomical" silicate grains at r = 35-110 AU, corresponding to the sub-mm ring plus an 
exterior halo, (3) a narrow belt at ~ 3 AU (T ~ 120 K) of small (a ~ 3 /jm) silicate grains, 
and (4) a narrow belt at ~ 20 AU (T ~ 55 K) of small grains (a ~ 8 fim) of undetermined, 
but possibly silicate, composition. Specific properties of the model components are presented 
in Table 2. 

Figure 7 displays the SED of the complete model and the separate disk components 
compared with the photometric and spectrophotometric data (model photosphere SED sub- 
tracted). Spitzer mid- and far-IR data especially reveal a complicated SED shape that 
strongly constrains the temperatures, locations, and grain sizes of warm unresolved mate- 
rial. The model is not unique but was built from the fewest components with the simplest 
assumptions that produced a good match to all the available data. Many alternate models 
were tested, resulting in confidence that: (a) 70 and 160 /zm emission does not extend be- 
yond 110 AU; (b) 350 /zm emission does not extend beyond 90 AU; (c) 350 /xm emission does 
not extend inside 35 AU; (d) the dominant emission in the sub-mm ring is not from silicate 
particles; (e) 70 and 160 /im emission within 35 AU must include non-emitting gaps, with 
emission restricted to annular zones; and, (f) grains in the innermost warm belt at r ~ 3 
AU have silicate composition. 

The model rises less steeply than the IRS observations at 10-18 /j,m. The excess flux 
in the model is less than 3% of the emission from the system in that range. This may be 
related to uncertainty in the photospheric model subtraction at short wavelengths where the 
IR excess is barely significant, and may also indicate the presence of crystalline silicates, 
with a sharp 20 /im spectral feature, combined with the amorphous silicates assumed in the 
model. 

The top panel of Figure 8 displays the model radial profile of perpendicular optical depth 
t_i_, i.e. fractional (unitless) surface density (m 2 of grain cross-sectional area per m 2 of disk 
surface area; e.g. Backman 2004) and in the bottom panel the corresponding mass surface 
density inferred from the optical depth and model grain properties. The slopes of the two 
functions across the sub-mm ring and halo are not equal because the model includes smooth 
variation in the particle size distribution between the fiducial radii, so the ratio of particle 
cross section (and thus optical depth) to mass surface density varies slightly with radius. 
The optical depths and mass densities of the two inner rings in Figure 8 are calculated by 
assuming Ar values of 1 AU at 3 AU and 2 AU at 20 AU; they would scale inversely with 
Ar. 

The total inferred mass in detected grains, i.e. with sizes from a few microns to a few 
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hundred microns, carrying most of the radiating surface area, is roughly 5 x 10~ 3 M e ~ 0.4 
MLuna, 90% of which is in the large ("icy") grain component of the sub-mm ring. The total 
bolometric fractional luminosity (L d /L*) of all the disk components is 1.1 xlO -4 . 

Figure 9 is a 2-D representation of the e Eri debris disk system in linear spatial scale, 
showing the sub-mm ring, halo, and inner warm belts. The position of the innermost warm 
belt is compared in the inset figure with one orbit solution for the candidate radial velocity 
planet (Section 5.4). 

5. Discussion 

5.1. Some Implications of the Disk Model 

The derived surface density profile for the small grain component (Figure 8) has a 
nearly constant value from 35 to 110 AU, consistent with dynamics controlled by Poynting- 
Robertson radiation drag or its corpuscular wind analog (Section 5.3). In contrast, the large 
grain component has surface density rising from 35 to 90 AU. The larger grains are less 
subject to drag forces, so their radial distribution may reflect that of the underlying parent 
body population. 

The void interior to the sub-mm ring is inferred to be partly filled by two additional 
belts of particles with mean sizes in both cases smaller than in the sub-mm ring. The large 
grain population in the ring appears to be absent in the void, even though 100 /zm-size ice 
grains could persist against sublimation for as long as e Eri's age, 850 Myr, at distances 
as small as 3 AU from the star (see equations 16 and 17 in Backman & Paresce 1993 for 
expressions of ice sublimation time scales). Moreover, the mid-IR SED, plus the 70 /im 
surface brightness profile in particular, indicate that there are gaps in the distribution of 
material from just beyond 3 AU to almost 20 AU, and from just beyond 20 AU to 35 AU. 
A mechanism is needed to prevent significant amounts of both small and large grains from 
drifting via P-R drag from the sub-mm ring into the inner void (Section 5.4). 

The total mass of detected grains in e Eri's innermost belt is about 0.1 x the mass 
estimated for the warm belt detected around the 230 Myr-old A star ( Lep (Chen & Jura 
2001). The size and structure of the e Eri system, with two warm belts in a central void 
surrounded by a cold disk, and silicate dust identified in the innermost ring, bear striking 
resemblance to the circumstellar material in the much younger HD 113766 system (10-16 
Myr, F3/F5 binary) (Lisse et al. 2008). 

An upper limit to the ratio of ice versus silicate grain mass for the sub-mm ring region 
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in which they co-exist is about 16. This ratio is intermediate between values for solar 
composition (~ 50) and for comets (~ 2-5, e.g. Altenhoff et al. 2002). 

5.2. Comparison with Prior Studies of e Eri 

5.2.1. Sub-mm Ring Features 

As noted above, the sub-mm ring at 350 /im (Figure 4) has the same inner and outer 
diameters observed in previously published 850 and 450 /iin maps (Greaves et al. 1998; 
2005), but 350 ring surface brightness features do not in general align with 450 or 850 
/im features. The (JCMT / SCUBA) 450 and 850 /im data were obtained via dithered 
("jiggle"-) mapping, whereas the (CSO / SHARC II) 350 /im map was made via continuous 
spiral scanning. Real features should be robust against such differences in data acquisition 
and reduction procedures. 

The strongest 450 and 850 /im feature, about 1.5 x brighter than the ring average, is 
in the south-east sector. A south-east feature appears in some form in all the maps and 
was suggested to be a real part of the e Eri system by Greaves et al. (2005) based on its 
motion over the 1997-2002 time span. Poulton, Greaves & Cameron (2006) made a careful 
statistical analysis of this issue. Only a weak feature is seen at that location in the 350 
/im image. Other fairly prominent 850 /im ring features, especially in the south-west sector, 
were considered by Greaves et al. to be background sources although they have no clear 
counterparts at 450 /tm. 

The strongest 350 /iin feature is located in the north-north-east section of the ring, 
with a surface brightness of about 0.24 mJy per square arcsec. Its location corresponds 
approximately to a weak 450 /im feature and a moderately strong 850 /itn feature that 
could also be part of the e Eri system, considering proper motion plus orbital motion. Recall 
(Section 3.1.3) that the 350 /tm data, when split into equal parts and composed into separate 
maps, have about the same level of disagreement regarding azimuthal features as between 
maps at different wavelengths, and between the 850 /iin maps for 1997-2002 versus 2000- 
2002 (Greaves et al. 2005). As discussed in Section 3.1.3 and displayed in Figure 2, no surface 
brightness features are seen in the ring at 70 and 160 /im. 

Real density features need not have the same surface brightness ratios versus wavelength 
as the main ring. For example, smaller grains trapped in planetary resonances can have 
increased libration widths, smoother distributions, and show less conspicuous clumping than 
large grains, so that some features prominent at long wavelengths may be less apparent at 
short wavelengths (Wyatt 2006). Azimuthal features with surface brightness that does not 
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vary progressively with wavelength, however, are difficult to reconcile with physical models. 

The disagreement between details of the different sub-mm maps illustrates the difficulty 
of working at wavelengths with large variable sky opacity and emission. These comparisons 
imply that the Spitzer far-IR images and CSO 350 /xm maps cannot confirm the ring surface 
brightness features reported at 450 or 850 fxm. 

5.2.2. Optical & Near-IR Upper Limits 
The brightness of scattered light from circumstellar material can be estimated from: 

Rx, r = (T ± /47r)(dl c /rl v )(P 0( 7 ) / cos(t) ) (2) 

in which R is the ratio, at a given wavelength and position, of circumstellar surface brightness 
per square arcsec to the star's brightness; t±_ is the perpendicular optical depth of the 
circumstellar material; d pc is the system distance from the observer; tau is the distance of 
the material from the star; P is the geometric albedo; 0(7) is the phase function value for 
the star-dust-observer angle 7; and cos(i) is the inclination of the disk to the plane of the 
sky, % = representing face-on presentation. A value of cos(i) ~ 1 can be used for the nearly 
face-on e Eri disk. A value 0(7) ~ 0.5 can be assumed for particles much larger than the 
wavelength of interest in an approximately face-on disk. 

Proffitt et al. (2004) used the STIS spectrographic camera on HST to search for scattered 
light from the sub-mm ring region. They set an upper limit of msTis <^ +25 for circumstellar 
light at r ~ 55 AU (17") in a wide bandpass with A e // ~ 0.6 /im. Assuming an albedo of 
P = 0.1 for the model large grains, which have the highest optical depth in the sub-mm 
ring region (Figure 6a), yields a predicted scattered light surface brightness at 55 AU almost 
exactly equal to the HST/STIS limit. Assuming instead an albedo of 0.9, possible for fresh 
pure ice surfaces, would yield an optical surface brightness almost 2.5 mags brighter than 
the limit. In comparison, the model of the 2nd warm belt (r ~ 20 AU ~ 6") predicts, an 
optical surface brightness of my ~ +21.7 per square arcsec assuming P = 0.1 and radial 
width Ar = 2 AU, brightness scaling inversely with Ar. 

Di Folco et al. (2007) report near-IR (K' band, 2.1 jum) interferometric observations of e 
Eri and r Ceti with the CHARA array. They resolved the stellar diameters in both cases and 
detected circumstellar emission near r Ceti but set an upper limit of 0.6% of the photosphere 
for emission within 1" (3 AU) of e Eri. The present model innermost belt predicts a total 
scattered light brightness of approximately 3% of the Di Folco et al. upper limit if the grain 



-24- 



near-IR albedo is assumed to be 0.1. The direct thermal emission at K' band from that belt, 
in fact from the entire disk, would be many orders of magnitude fainter than the scattered 
light brightness. 

5.2.3. SED Models 

Dent et al. (2000), Li, Lunine & Bendo (2003) (hereafter LLB03), and Sheret, Dent 
& Wyatt (2004) (hereafter SDW04) combined pre-Spitzer broadband photometry and sub- 
mm map data to produce models of the e Eri disk with specific grain spatial distributions, 
size distributions and compositions. Given the limited data available to constrain them, 
those models are consistent with the model presented in Section 4 of this paper based on 
higher-precision observations over a much broader wavelength range. 

Dent et al. (2000) assumed non-porous grains, with a simplified emissivity law deter- 
mined only by particle size without mineralogical features. Their model requires material 
as close as 10 AU from the star with ~ 10% of the surface density in the main ring. LLB03 
and SDW04 calculated models assuming real grain mineralogies, including a range of porosi- 
ties (fraction of silicate grain volume occupied by vacuum or ice) ranging from p = to 
0.9. SDW04 included their own sub-mm photometric data regarding e Eri to supplement 
published maps. LLB03's best-fit model involved highly porous (p = 0.9) grains, whereas 
SDW04's suite of models favored non-porous (solid) grains. The present model assumes 
non-porous grains and agrees qualitatively with the SDW04 result that the grain size distri- 
bution cannot extend down to 1 /im without producing too much short-A emission. SDW04 
remarked that the lack of /zm-scale grains is a puzzle given the inability of e Eri's radia- 
tion field to eject grains of that size. Radiation pressure and corpuscular wind effects are 
discussed further in Section 5.3. 

The total circumstellar grain mass derived by LLB03 was ~7x 10~ 3 M e ~ 0.6 M^ una , 
compared with 0.1 M.L una in SDW04's model and 0.4 M^ una in the current model. The 
differences between these mass estimates result from differing assumptions about the particle 
properties and size ranges, since each group is matching approximately the same system 
fluxes, requiring similar total particle radiating areas. 
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5.3. Particle Dynamics 



5.3.1. Grain Ejection via Wind versus Radiation Pressure 



Small grains are subject to radiation pressure (e.g. Burns et al. 1979) and corpuscular 
winds (e.g. Plavchan et al. 2005) that tend to remove them from their source region, either 
rapidly outward by direct ejection for particles smaller than the "blowout" size, or, for larger 
particles, slowly inward toward encounters with planets or sublimation near the star. Grains 
also collide destructively. Competition between rates of destructive processes versus rate of 
grain resupply from larger parent bodies determines a debris disk's equilbrium density and 
also whether there is net mass flow outward or inward from the grain source (e.g. Wyatt 2005; 
Meyer et al. 2007). These processes will be considered in this and the following subsection in 
view of the e Eri system's characteristics. Porous and non-spherical grains would be affected 
more strongly by radiation and wind than the spherical and non-porous grains assumed here, 
so destruction and drift time scales calculated below should be considered upper limits. 

The direct (radial) radiation pressure (force per area) opposing gravitational attraction 
between star and grain is: 



The wind mass loss rate for e Eri has been determined to be ~ 3 — 4 x 10 10 kg s _1 (Stevens 
2005; Wood et al. 2005), about 20-30 x the solar value. Assuming conservatively that young, 
chromospherically active e Eri has the same wind terminal speed of ~ 400 km/sec as the 
sun's, radial wind pressure is negligible in the e Eri system, less than 5% as strong as radiation 
pressure. The star's mass/luminosity ratio is low enough that solid silicate-density particles 
should be stable against direct ejection down to sizes of a few x 0.1 /xm. 

This analysis indicates that: (1) the model grain size lower limits of a few microns in 
the sub-mm ring and inner unresolved belts are not defined by radiation pressure or neutral 
wind "blowout", and (2) radiation and wind pressures in this system are too weak to move 
substantial amounts of material outward, for example from the sub-mm ring into the halo. 



P r (radiation) = L*/47rr 2 c 



(3) 



whereas the direct wind pressure is: 




(4) 
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5.3.2. Poynting-Robertson and Wind Azimuthal Drag versus Collisions 

The following two equations compare the azimuthal pressure Pq on orbiting grains re- 
spectively due to radiation (P-R drag) and a neutral wind. 

Pg{radiation) = v or bL*/4:irr 2 c 2 (5) 

Peiwind) = v orb (dMjdt)/Arrr 2 (6) 

Using the values given above for wind mass and speed, the radial drift rate for particles 
around e Eri is enhanced by a factor of about 20 relative to the purely radiation-driven P-R 
drag rate (numerically similar, by coincidence, to the e Eri wind mass loss rate relative to 
solar). 

The timescale for grain-grain mutual collisions is: 

tcoll ~ Porb/ST^ (7) 

in which P or b is the orbit period and tj_ is the perpendicular optical depth. Cross-velocities 
as low as a few hundred meters per second are sufficient to destroy small silicate and ice 
grains. The solar system's asteroid and Kuiper belts both have inclination and eccentricity 
ranges more than large enough to result in routinely shattering collisions, which may also 
be assumed to pertain to the e Eri disk. 

At the inner edge of the sub-mm ring, r = 35 AU, where P-R drag is most competitive 
with collisions, the collision time scale for the large (a ~ 135 /im) grains is t co u ~ 2 x 10 5 yrs, 
whereas the wind-drag time scale for inward drift by Ar = 1 AU is ~ 8 x 10 5 yrs. Thus, even 
when enhanced by wind, virtually all the large icy grains should be destroyed by collisions 
in situ and not be able to drift inward to fill the central void. For the small grains (a ~ 15 
/im), on the other hand, the collision destruction time scale (including cross-sectional area 
of both small and large grains as targets) is ~ 1.5 x 10 5 yrs, approximately the same as 
the wind-drag time scale for these particles to drift 1 AU. Given that the collision process 
is statistical, e.g. ~ 1% of grains can be expected to survive 5 t co u, some fraction of the 
small grain population could travel inward from the sub-mm ring into the central void before 
colliding with other grains. 

At the inner and outer edges of the halo, r = 90 and 110 AU respectively, the collision 
time scale for the small grains is also comparable to the Ar = 1 AU wind-drag drift time. 
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Thus, if grains are produced in the halo region they would tend to stay in place or flow 
inward toward the sub-mm ring, rather than outward. 



Is there a mechanism other than radiation or wind pressure that could move material 
from the sub-mm ring to the outer halo, or does the presence of significant material in the 
halo require a particle source (presumably collisions of planetesimals) at r > 90 AU? In fact, 
scattering of material into new orbits after collisions could suffice. 

The vis viva equation expresses the orbit properties resulting from velocity other than 
circular orbit velocity at a given location, 



in which v is the speed at radius r, and a or b is the orbit semi-major axis. 

Under what conditions would grain collisions in the main, heavily populated ring be 
able to supply material to the outer halo? If grains or parent bodies collide at r with average 
cross-speeds of £ x v C i rc , corresponding to relative orbit inclinations of order 5i ~ ±£ rad 
and/or eccentricity ranges of order 8e ~ £/2, significant numbers of fragments will leave 
the collision with 5v > £v relative to the motion of the colliding bodies' center-of-mass. 
Collisions at r ~ 90 AU with £ = 0.05 would send some fragments into new orbits reaching 
to 110 AU, the outer edge of the halo. Therefore the extended halo observed at 70 and 
160 jum could be the result of collisions among planetesimals and fragments moving with 
modest departure from circular planar motion in the sub-mm ring. If collsions are in fact 
the mechanism spreading material from the ring into the halo, the lack of material beyond 
110 AU could indicate an upper limit on inclinations and eccentricities in the ring. 



Radiation plus wind pressure should not be able to eject solid silicate-density particles 
larger than sub-micron size from the e Eri system, but the disk model has a deficit of particles 
smaller than a ~ 6 /im in the sub-mm ring. Consideration of particle destruction time scales 
suggests that a time span of order 10 5 years may have elapsed since the last significant 
injection of small particles. The small grain population has a spatial distribution indicating 
overall control by azimuthal drag. Some small grains near the inner edge of the sub-mm ring 



5.3.3. Collisional Scattering 




(8) 



5.3.4- Particle Dynamics Summary 
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should be able to drift into the central void faster than their collision destruction time scale. 
The fact that this does not appear to be happening may indicate the presence of a barrier 
planet preventing inward drift (Section 5.4). 

The large grains in the sub-mm ring have collision time scales much shorter than drift 
time scales, so they are seen near their locations of creation, and their radial mass density 
profile may track that of the planetesimal parent bodies. 

5.4. Planets Around e Eri 

Radial velocity measurements of e Eri have suggested the presence of a planetary com- 
panion, hereafter designated Planet A, with M sini = 0.86 M Jup and an orbit period of about 
6.9 years (Hatzes et al. 2000 and references therein). The corresponding orbital semi-major 
axis is 3.4 AU, near the location inferred for the innermost warm dust belt. 

The eccentricity of 0.70 ± 0.04 reported for the 3.4 AU planet by Benedict et al. (2006) 
corresponds to periastron and apoastron distances of 1.0 and 5.8 AU, respectively. A giant 
planet with this orbit would quickly clear that region not only of dust particles but also the 
parent planetesimal belt needed to resupply them, inconsistent with the dust distribution 
implied by the Spitzer observations. The disk's 70, 160 and 350 /iin aspect ratios and position 
angles reported here do not correspond to values of those quantities determined from the 850 
/iin map, calling into question the system inclination value used to combine the astrometric 
and radial velocity data. Butler et al. (2006) report a different orbital eccentricity of 0.25 ± 
0.23, which could allow the planet to orbit entirely outside the innermost dust belt. 

As discussed in Section 5.3, some of the small grains in the sub-mm ring should be 
able to drift into the central void and begin to fill it on timescales of 10 6 to 10 7 years. A 
mechanism is needed to prevent grains from moving inward past 35 AU. Liou & Zook (1999) 
and Moro-Martin et al. (2005) modeled the ability of massive planets to define the inner 
edges of Kuiper belt-like dust rings. A hypothetical planet orbiting at the edge of the main e 
Eri sub-mm ring is here designated Planet C. The gravitational influence of a planet near or 
in the ring has also been proposed (Liou & Zook 1999; Quillen & Thorndike 2002; Greaves 
et al. 2005) as the cause of ring density enhancements corresponding to surface brightness 
features reported in sub-mm imaging that may or may not be real (Section 5.2.1). 

Confinement of a narrow debris belt at r ~ 20 AU could indicate the presence of a third 
planet, here designated Planet B. Deller & Maddison (2005) calculated dynamical models of 
the e Eri system including effects of planetary masses, radiation pressure, and stellar wind, 
and found that a planet orbiting at the inner edge of the sub-mm ring would not by itself 
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be able to define the ring's edge. Rather, the effect of an additional Jovian-mass planet 
with semi-major axis a ~ 20 AU is needed. Greaves et al. (2005) proposed that the mass 
controlling some of the possible sub-mm ring features may be located well inside the ring, 
at r ~ 24 AU, based on the orbital period suggested by apparent motion of those features 
over 7 years. 

Summarizing all the observations and models bearing on the presence of planets around 
e Eri, there may be indications of three different planets in the system. Planet A: the 
long-suspected but still unconfirmed Jovian-mass planet with a or b = 3.4 AU that may be 
associated with the innermost warm debris belt detected by Spitzer. Planet B: also perhaps 
Jovian-mass, associated with the second warm debris belt at r ~ 20 AU inferred in this 
paper, also possibly helping define the sub-mm ring inner edge via 2:1 resonance, keeping 
the zone between 20 and 35 AU clear. (Note that a 1 M Jup planet with age 0.85 Gyr located 
6" from e Eri would not have been seen by Spitzer but could be detected by JWST.) Planet 
C: located at r ~ 35 AU, with mass less than a few x 0.1 Mj up according to dynamical 
models, preventing small grains from drifting inward past the inner edge of the sub-mm 
ring. 



5.5. Simple Collisional-Evolutionary Model 

A dynamical model of the e Eri sub-mm ring was calculated assuming equilibrium 
between collisional production of grains and removal by P-R plus corpuscular wind drag, 
based on previous simple dynamical models of the solar system's Kuiper belt (Backman, 
Dasgupta & Stencel 1995) and solar-type circumstellar debris disks observed by Spitzer 
(Meyer et al. 2007). This model indicates that the sub-mm ring can be maintained by 
collisions of 11 M ffi of parent bodies with assumed diameters of 10 km and material densities 
of 1000 kg m~ 3 . This agrees with estimates made by Sheret, Dent & Wyatt (2004) and 
Greaves et al. (2005) using similar assumptions. Thus, the e Eri system as observed can be 
in quasi-equilibrium. 

Spitzer observations of Vega revealed a far-infrared disk extending across the position 
of a sub-mm ring into an outer halo (Su et al. 2005), similar to the e Eri system. However, 
the Vega disk's properties require recent injection of circumstellar material to explain the 
abundance of short-lived particles. The main difference between the two systems is that e 
Eri's luminosity (including wind) is low relative to Vega's, so there is no direct ejection of 
grains, P-R plus wind azimuthal drag rates are lower, and equilibrium resupply by continuous 
minor collisions can be effective. 
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A simple evolutionary model that estimates loss of belt mass via destruction of plan- 
etesimal parent bodies, with subsequent removal of fragments by the dynamical processes 
listed above, was described and applied to the solar system's Kuiper belt, asteroid belt, and 
extrasolar analogs by Meyer et al. (2007). The same model, tuned to the properties of the e 
Eri system, predicts that for an initial mass of planetesimal parent bodies in the range 15-30 
M ffi , similar to the inferred original mass of the Kuiper belt (e.g. Stern & Colwell 1997), 
about 11 M e would remain at a system age of 0.85 Gyr. In other words, e Eri's current 
debris disk is similar in total mass to what the solar system's Kuiper belt mass would have 
been at the age of e Eri had there not been a prior Late Heavy Bombardment (LHB) clearing 
event (e.g. Levison et al. 2007). 

Collisional evolution is dissipative, eventually converging on a planetesimal population 
that has a time scale for further collisions proportional to the system age, nearly independent 
of the initial conditions, if there are no external effects. However, to transform the young 
solar system's Kuiper belt or e Eri's debris disk into a low-mass system like the present-day 
Kuiper belt, with at most a few x 0.1 M ffi at age 4.5 Gyr (Luu & Jewitt 2002), requires 
additional mass reduction processes beyond quasi-steady state collisional grinding and debris 
removal (Meyer et al. 2007 and references therein). For example, in our solar system, the 
LHB event, possibly driven by a late episode of planet migration, may have destroyed most of 
the Kuiper belt's mass relatively sudddenly, and similar events may happen in other systems 
with ages of a Gyr or more (Levison et al. 2007 and references therein). Alternatively, 
a process that is continually at work in our solar system, planetary pertubations pulling 
planetesimals out of the Kuiper belt, could, with the help of collisional grinding, transform e 
Eri's 11 M e disk at t — 0.85 Gyr into the solar system's 0.1 M e Kuiper belt at 4.5 Gyr if the 
perturbative removal process supplementing collisions were to operate with an exponential 
time scale of about 1.5 Gyr. In contrast, r Ceti's debris disk has an inferred mass of 1.2 M e 
at an age of 8 Gyr, so it may not have had a major disruptive episode (Greaves et al. 2004). 



6. Summary 

The planetary debris disk around e Eridani was observed using Spitzer's IRAC and 
MIPS cameras, Spitzer's IRS and MIPS SED-mode spectrophotometers, and the SHARC II 
bolometer array on the CSO sub-millimeter telescope. 

The slight 25 /iin excess detected by IRAS is confirmed photometrically by the Spitzer 
IRS and MIPS data but not resolved in the 24 /zm image. The CSO 350 /mi map confirms 
the sub-mm ring at r ~ 11-28" (35-90 AU), previously observed at 450 and 850 jum, about 
2x the width of our solar system's Kuiper belt. Analyses of the 70 and 160 jum images 
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reveal emission extending from within a few arcsec of the star outward to r ~ 34" (110 
AU), both inside and outside the sub-mm ring. The Spitzer 70 and 160 /zm images do not 
have enhanced surface brightness at the position of the ring. The 350 /j,m surface brightness 
varies noticeably with azimuth around the ring but the locations of the brightest features 
do not correspond between 350 /zm data sets reduced separately, or with with the brightest 
features in the longer-wavelength maps, calling into question the reality of the sub-mm ring 
azimuthal asymmetries. 

The over-all system spectral energy distribution is complex, with significant excess rela- 
tive to the photosphere beginning at A ~ 15 /zm. The excess rises rapidly to 20 /im, flattens 
from 20 to 30 /im, rises rapidly again from 30 /im to a peak between 70 and 100 /im, and 
then decreases toward longer wavelengths. 

A model of the disk, constrained by all the available data, combines two populations of 
material in the resolved disk: (1) large (radii a ~ 100-200 /im) solid grains in the sub-mm 
ring, which have high emissivity at sub-mm wavelengths but low emissivity at 70 and 160 
/im, consistent with H 2 ice composition, plus (2) smaller (a ~ 6-23 /im) solid grains in 
both the sub-mm ring and its exterior "halo" (r = 25-34", 90-110 AU) which have high 
emissivity in the far-IR and low emissivity at sub-mm wavelengths, consistent with silicate 
composition. The mid- and far-IR spectrophotometry and radial source profiles imply two 
additional narrow belts of material, one at r ~ 3 AU (1") with a weak 20 /im silicate emission 
feature, and the other at r ~ 20 AU (6") of undetermined, but possibly silicate, composition. 

The grains in the halo may be spread there from the ring region via the effect of collisions. 
Definition of the inner edge of the sub-mm ring, and confinement of the two rings of material 
closer to the star, probably require gravitational influences of three planets. A dynamical 
model of the disk indicates that the observed populations of grains in the sub-mm ring and 
outer halo can be supplied in equilibrium by collisions of about 11 M e of km-scale parent 
bodies located in the sub-mm ring. 

This paper is based in part on observations made with the Spitzer Space Telescope that is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under NASA 
contract 1407, and in part on observations from the Caltech Submillimeter Observatory 
(CSO). Support for this work was provided by NASA through awards issued by JPL/Caltech, 
including subcontracts 1278243 and 1255094 respectively to the SETI Institute and the 
University of Arizona, and also by NASA grant NNG05GI81G. Analyses in this paper made 
use of the SIMBAD database operated at CDS, Strasbourg France. 

We are grateful to Jane Morrison and Myra Blaylock, University of Arizona, for help 
with MIPS image data processing, and Paul Smith, University of Arizona, for guidance in 
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Table 1. e Eri Broadband Photometry 



A [/xm] 


Total F v 


[Jy] 


Star F v [Jy] 


Excess F v [Jy] 


Notes 


3.550 


64.320 ± 


1.140 


64.320 




IRAC 


4.493 


39.800 ± 


0.700 


39.800 




IRAC 


5.731 


26.110 ± 


0.470 


26.110 




IRAC 


7.872 


14.910 ± 


0.260 


14.910 




IRAC 


12.000 


6.690 ± 


0.540 


6.570 


0.130 


Color-corrected IRAS 


23.680 


2.042 ± 


0.040 


1.726 


0.315 


MIPS 


25.000 


1.900 ± 


0.280 


1.550 


0.350 


Color-corrected IRAS 


60.000 


1.750 ± 


0.250 


0.270 


1.480 


Color-corrected IRAS 


71.420 


1.688 ± 


0.207 


0.189 


1.499 


MIPS 


100.000 


1.990 ± 


0.300 


0.100 


1.890 


Color-corrected IRAS 


156.000 


0.970 ± 


0.186 


0.039 


0.931 


MIPS 


200.000 


1.880 ± 


0.280 


0.020 


1.860 


ISO 1 


350.000 


0.366 ± 


0.050 


0.008 


0.358 


CSO / SHARC II 


450.000 


0.250 ± 


0.020 


0.005 


0.245 


JCMT / SCUBA 2 


850.000 


0.037 ± 


0.003 


0.001 


0.036 


JCMT / SCUBA 2 



References. — (1) Walker & Heinrichsen (2000); (2) Greaves et al. (2005) 
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Table 2. Model Components 



component r [AU] 


Mj 




] 


a 


a [pan] 


X 




/ 




Wl 


3 


1.8 


X 


10" 


-7 




3.0 




3.3 


X 


10" 


-5 


W2 


20 


2.0 


X 


10" 


-5 




8.0 




3.4 


X 


10" 


-5 


RS 


35-90 


2.0 


X 


10" 


-4 


+0.01 


6.0-23 


-3.5 


3.0 


X 


io- 


-5 


RL 


35-90 


4.2 


X 


io- 


-3 


+ 1.05 


100-200 


-3.5 


4.4 


X 


io- 


-6 


HS 


90-110 


2.5 


X 


10" 


-4 


+0.15 


15-23 


-3.5 


4.8 


X 


io- 


-6 



Note. — Columns: (1) model component: Wl = warm belt 1, W2 = warm belt 2, RS = 
sub-mm ring, small grains; RL = sub-mm ring, large grains; HS = halo, small grains; (2) 
location; (3) total mass; (4) mass surface density exponent, assumed to be zero for the Wl 
and W2 components, fitted to data for the other components; (5) grain radius; (6) assumed 
grain size distribution exponent; (7) fractional luminosity, Ld/L* 
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Fig. 1. — Images of e Eri in the three Spitzer MIPS channels. N is up and E to the left in 
all the panels; 10" = 32 AU at e Eri's distance. A cross marks the stellar position at the 
epoch of each observation, taking into account e Eri's proper motion of ~ 1" per year. Dark 
is positive. Circles indicate FWHM beam sizes. The left column shows 24 /j,m images 
in log stretch. These have been reconstructed from multiple dithered exposures and are 
oversampled by a factor of four. Panel (a) is the direct image, dominated by the stellar PSF. 
Panel (d) shows this image after subtracting a PSF reference star scaled to the photospheric 
flux density. The remaining excess (18% of the total 24 /im emission) appears largely as an 
unresolved point source. The middle column shows 70 /im images in linear stretch. Panel 

(b) is the default pixel-scale image used for photometric measurements, and panel (e) shows 
the fine pixel-scale image. The right column shows 160 /im images in linear stretch. Panel 

(c) is the direct image, and panel (/) shows this image after subtraction of the ghost image 
arising from the spectral leak. 
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Fig. 2. — (a) difference between the Spitzer MIPS fine-scale 70 /xm image and its azimuthal 
average, (b) the same for the MIPS 160 fxm image. The vertical bars are due to incomplete 
correction of column artifacts. The remaining residuals are consistent with noise from the 
subtraction, indicating absence of surface brightness features in the sub-mm ring within the 
sensitivity limits. 




Fig. 3. — MIPS SED-mode measurements of e Eri. (a) (solid line) FWHM of the e Eri disk 
emission versus wavelength, extracted by Gaussian fitting along the slit at PA 85°, versus 
(dashed line) the same for Canopus, an unresolved calibration star with no IR excess, e 
Eri is clearly resolved across the wavelength range of the MIPS SED measurements, (b) 
Star+disk spectrum windowed by the 20"-wide slit, without aperture correction. The dotted 
lines indicated the range of photometric uncertainty based on uncertainty in Gaussian fits 
across the slit. Subsequent model-dependent flux normalization of these data is discussed in 
Section 4.4.2. 
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Fig. 4.- Caltech Submillimeter Observatory (CSO) SHARC II 350 fim map. Dark is 
positive. The circle's diameter indicates 12.7" resolution to which the map was smoothed, 
equivalent to 41 AU at the distance of e Eri. A scaled PSF representing the relatively small 
350 /xm brightness of the central star has been subtracted. 
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Fig. 5. — Merged broadband SED of the e Eri system. Numerical values are presented in 
Table I. The thin solid curve is a Kurucz photosphere for T = 5000 K, log g = 4.5, fit to 
the four Spitzer IRAC photometric points at 3.6 to 7.9 /im. The thick curve is the combined 
Spitzer IRS SL + SH + LH spectrum rescaled to have zero average excess at 5-12 /im. The 
IRAS 100 /im point and the ISO 200 /im point are considered upper limits due to background 
contamination. The inset panel shows the mid-infrared segment of the SED from Spitzer IRS 
(thin solid line), after subtraction of the photosphere model displayed in the main figure. 
The thick line is the same spectrum smoothed with a Ill-point (~ 3.3 /im) boxcar. 
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Fig. 6. — Radial profiles of images at 24, 70, 160 jttm (Spitzer MIPS) and the 350 jum (CSO) 
map. In each wavelength panel the dashed line with error bars represents the azimuthal 
average of the image data after subtraction of the stellar photosphere PSF. Observational 
error bars are estimated as the rms of the azimuthal average. The thick solid line is the disk 
model after convolution with the relevant instrumental PSF, fitted to the observations. The 
discontinuous thin solid line is the model intrinsic surface brightness before PSF convolution 
(Section 4.2). 
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Fig. 7. — Observed SED of the e Eri disk, after subtraction of the stellar photosphere, com- 
pared with the model SED. Individual photometric points include MIPS SED flux densities 
at 55 and 90 /xm, iteratively aperture-corrected using model images, and rescaled according 
to the measured MIPS 70 /im image total flux. The thick solid line is the IRS SL + SH + 
LH combined spectrum, rescaled to have zero average excess for A = 5 — 12/im. The thin 
solid line is the total model flux (sum of all four dust components). The two dashed lines are 
the contributions of the two unresolved inner belts. The dot-dashed line is the contribution 
of the small (a ~ 15 /im) silicate grains in the sub-mm ring and halo. The dotted line is the 
contribution of the large (a ~ 135 /im) ice grains in the sub-mm ring. 
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Fig. 8. — (a) Perpendicular optical depth t± versus position for the five model disk com- 
ponents. The two unresolved inner belts (dashed lines), outer "icy" ring (dotted line), and 
ring+halo " silicate" dust populations (dot-dashed line) are shown, (b) Mass surface density 
versus radius corresponding to the optical depth profiles in panel (a), assuming grain com- 
positions and material densities described in the text. The difference in slopes between the 
optical depth and mass surface density profiles is due to radial variation in model particle 
properties. 
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Fig. 9. — Representation of the e Eri debris disk model components. The small-scale dotted 
ellipse is one solution for the orbit of the suggested radial velocity planet that appears to be 
inconsistent with the innermost warm debris belt's position. 



